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Abstract 

The entanglement in one-dimensional random XY spin systems where the impurities of exchange cou- 
plings and the external magnetic fields are considered as random variables is investigated by solving the 
different spin-spin correlation functions and the average magnetization per spin. The entanglement dynam- 
ics near particular locations of the system is also studied when the exchange couplings (or the external 
magnetic fields) satisfy three different distributions(the Gaussian distribution, double-Gaussian distribu- 
tion, and bimodal distribution). We find that the entanglement can be controlled by varying the strength of 
external magnetic field and the different distributions of impurities. Moreover, the entanglement of some 
nearest-neighboring qubits can be increased for certain parameter values of the three different distributions. 
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I. INTRODUCTION 



Entanglement not only has the interesting properties of quantum mechanics but also is very 
important in the quantum information processing (QIP), such as quantum teleportationM, dense 
coding PI, quantum secret sharing ' 3 1, quantum computation W and some cryptographic protocols' 5 '. 
In order to realize the quantum information process, a great effort has been made to study and 
characterize the entanglement in cavity QED' 6 7 ' and solid state systems. A typical example is the 
spin chains that can describe interaction of qubits not only in solid physical systems but also in 
many other systems such as quantum dots' 8 ', electronic spins' 9 ', and optical lattices' 10 ]. Therefore, 
there have been conducted numerous studies on Ising model' 11 ' and all kinds of Heisenberg XY 
XXZXYZ models' 12 - 15 !. 

Impurities often exist in solid systems and play an important part in condensed matter physics. 
As a candidate of QIP, a solid system with impurity is also one of our important study objects. In 
the previous researches, the effect of impurity on the quantum entanglement has been studied in 
a three-spin system I 16 " 17 ' and a large spin systems under zero temperature' 18 '. However, in these 
studies, only single impurity has been studied. 

Recently, Huang et al' 19 20 ' have demonstrated that for a class of one-dimensional magnetic 
systems entanglement can be controlled and tuned by varying the anisotropy parameter in the XY 
Hamiltonian and by introducing impurities into the systems. However, in Ref.(19), only the impu- 
rity and the external magnetic fields in a Gaussian form are considered and the value of the width 
of the distribution is fixed. In Ref.(20), the strength of the impurity is located at two sites. For 
the pure case, Osterloh et al' 21 ' examined the entanglement between two spins at position i and j 
in the spin chains. Owing to its importance, in this paper we study the entanglement dynamics 
near particular locations of one-dimensional s — \ random XY spin system when the exchange 
couplings (or external magnetic fields) satisfy three different distributions(the Gaussian distribu- 
tion, double-Gaussian distribution, and bimodal distribution), to our knowledge, which have not 
been reported yet. The present study in simple examples can help us to understand the behaviour 
of the entanglement in one-dimensional random XY spin systems for the different distributions. 
More importantly, we will demonstrate that one can control or manipulate the entanglement in 
spin system with the help of the exchange couplings and the external magnetic fields. 
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II. SOLUTION OF THE XY MODEL 



We consider a physical Heisenberg XY model of N spin-| particles interacting with their near- 
est neighbours. In the presence of impurities, the one-dimensional Hamiltonian is given by I 19 l 



N 



1+7 x x 1-7 

H = y~ Ji ^+ ia i a i+i — y~ 

N N 



1=1 1=1 

where J^+i is the exchange interaction between sites i and i+1, hi is the strength of the external 
magnetic field on site i, a x,y,z are the Pauli matrices, 7 is the degree of anisotropy and N is the 
number of sites. The periodic boundary conditions satisfy cr x N+1 = of, <? y N+l = of, cr z N+1 — a i- 
Now we define the raising and lowing operators af ,a~ and introduce Fermi operators I 22 ' c+ and 
Cj, and they are expressed as follows: 

1 l ~ l 
4 = 2 {< + = 4 exp(i7r ^ ct Cj ) , (2) 

j'=i 

1 1-1 

j'=i 

so that, the Hamiltonian has the following form 

# = - Ji,i+i[(4c i+ i + 7c, + c+ J + h.c] 
i=i 

TV 

- 2 J>(c+ Ci --) (4) 

i=i 

In the present paper, the exchange interaction has the form J i>i+ i — J(l + a^j+i), where a intro- 
duces the impurity in the double-Gaussian form with peaks at with strength d and at ^f^- 
with strength ( 2 , 

N + 1 

= P x Ci exp{-e(? — )} 

+ (1 - P) x C 2 exp{-e(* - ^-*)} (5) 

The external magnetic field takes the form hi — h(l + Pi), where 

iV + 1 

# = P x & exp{-e(i — )} 

N - 1 

+ (1 - P) x &exp{-e(i _ — — )} (6) 
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When a = (3 = 0, the above reduces to a pure case; when P — 1, the above reduces to the case in 
Ref.(19). By introducing the dimensionless parameter A = J/ (2h), the symmetrical matrix A and 
the antisymmetrical B, the Hamiltonian becomes 



N 1 
H = E + ~(4Bi,jC+ + h.c)\ 



(7) 



The above Hamiltonian can be diagonalized by making linear transformation of the fermionic 
operators r\ k = J2iLi 9ki c i + h ki cf, r]£ = J2iLi 9kicf + h ki Ci, and then the Hamiltonian becomes 



JV 



H = ^ A kVtVk + const, 



(8) 



k=l 



and two coupled matrix equations satisfy 4> k (A — B) = Ak4>k,4>k(A + B) = A k (p k , where the 
components of the two column vectors <p ku tp ki are given by <p ki = g ki + h kh i/j ki = g ki - h ki . 
Finally, the ground state of the system \ip ) can be written as r/ k \tp ) = 0. 

III. SPIN-SPIN CORRELATION FUNCTIONS 

Before we dicuss the entanglement, we should have a brief review of spin-spin correlation 
functions. The spin-spin correlation functions for ground state and the average magnetization per 
spin are respectively defined as ^ 

S lm = \ fa Wf<\ ^0) , Sf m = \ fa |(7f (7* | ^0> , 

SL = \ fa Wf<\ ^o) , Mi = \ \o-tfa) . 



These correlation functions are given as expectation values of products of fermion operators. 
Using Wicks theorem f 23 ', these expressions can be rewritten as 



ox 1 

D lm — 4 



qy _ l 



^1,1+1 ^1,1+2 ■ ' ■ ^l,m 



y G m -lJ+l Gm-1,1+1 ■ ■ ■ Gto-1,to j 

Gi + ij Gi+ij+i ■ ■ ■ Gi + 

l,m— 1 1 

y G m ,i G m j + i ■ ■ ■ G m>m -i J 
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Sim — l{Gl,lGm,m ~ G m jG^ m ), Mf — ^G 



where G 



h3 



kj 



IV. ENTANGLEMENT OF NEAREST-NEIGHBOURING QUBITS 

In this part, we give the expression of the concurrence that quantifies the amount of the entan- 
glement between two qubits. For a system described by the density matrixp, the concurrence C 

is [24] 



C(p) = max(0, Ai - A 2 - A 3 - A 4 ) 



(9) 



Here Ai, A 2 , A 3 , A 4 are the eigenvalues (of them Ai is the largest) of the spin-flipped density 



operator R, which is defined by R = yfppyfp , where p = (a y <S> o- y )p*{a y <g> a y ), p denoting the 
complex conjugate of p with a y being the usual Pauli matrix. The values of concurrence C ranges 
from zero to one; when C — 0, the two qubits are in an unentangled state, when C = 1, the two 
qubits are in an maximally entangled state. 

Using the operator expansion for the density matrix and the symmetries of the Hamiltonian' 25 ], 
in the basis states , \][) , |J,f) , |||)}, p has the general form 

( Pi,i p 1A \ 

P2,2 P2,3 

P3,2 P3,3 

\p 4 ,l p 4 , 4 / 
with 

Aa = VPl,lP4,4 + |pi, 4 | , A 6 = y/p2,2p3,3 + 1^2,3 I , 
Ac = y/pTJp~44 - |P1,4| , A d = v /p 2 ,2P3,3 - 1^2,3 1 , 

We can express all the matrix elements in the density matrix in terms of different spin- spin 
correlation functions: 



Pi,i = \Mf + \M z m + Sf m + i, 
P2, 2 = \Ml - \M* m - Sf m + J, 
p^ = \M* m -\Mt-Sl m + \, 
P4,4 = -\Mf - + S* m + \, 
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P2,3 - Sf m + Sf m , 

„ Qx cU 

Pl,4 — O lm - O lm , 



V. RESULTS AND DISCUSSIONS 



In this paper, we focus our discussions on the transverse Ising model with 7 = 1. Our goal is to 
examine the dynamics of entanglement in the varying of the exchange couplings and the external 
magnetic fields. First, we examine the change of the entanglement for the nearest neighbouring 
concurrence C(i,i+l)for different values of the impurity as the parameter A varies. We consider 
two kinds of nearest neighbouring concurrences near particular locations of the system. Figure 
1 depicts the nearest neighbouring concurrence C(49,50) as a function of the reduced coupling 
constant A at different values of the impurity ( for different distributions with the system size N 
=101 and the anisotropy parameter 7=1. Figure 1(a) shows the change of concurrence C(49,50) 
as a function of different values A with p — 1, i.e the Gaussian distribution. We can see that the 
concurrence increases and arrives at a maximum close to the critical point A c , and it is close to 
zero above A c . As ( increases the concurrence tends to increase faster and the A m , where concur- 
rence approaches a maximum, shift to left very rapidly. This is consistent with the result in Ref. 
[19](Fig. 1). In Fig. 1(b), 1(c), and 1(d), we give the curves for the concurrence against the width 
of the double-Gaussian distribution. The two Gaussian distributions have equal probability with 
p = 0.5, and the central positions are at and ^-k Here one of the double-Gaussian distribu- 
tion is fixed with d = 0.5. We first investigate the situation when the width of the double-Gaussian 
distribution e is 0. 1 , A similar behaviour can be seen in Fig. 1 (b), only the changed width becomes 
narrow. As e increases, the concurrence increases slowly and the peak value decreases, which 
is shown in Fig. 1(c). As is well known, bimodal distribution is a particular case of the double- 
Gaussian distribution, that is to say, the double-Gaussian distribution is converted into bimodal 
distribution as e increases. In Fig. 1(d), e = 10. The numerical calculations show that concurrence 
decreases with the increase of ( 2 , which indicates that the behaviours are very different from the 
former cases. 

In Fig. 2, we show the results of the nearest neighbouring concurrence between the sites 49 and 
50, as a function of the parameter A for different strengths of the external magnetic field £. The 
effect of the external magnetic field £ in the Gaussian distribution is also shown in Fig. 2(a). How- 
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ever, different from the effect of the exchange couplings, the concurrence increases slowly and 
tends to move to infinity by increasing the value of the parameter £. This is also consistent with 
the result in Ref. [19](Fig. 1). A similar behaviour can be seen in Fig. 2(b) for the double-Gaussian 
distribution, however, with £ = 1 and £ = 10, it is interesting to find that the entanglement peak 
between the nearest neighbours increases to a value larger than that in Fig. 2(a). As £ increases, 
the concurrence increases rapidly below A c , while the concurrence increases slowly above A c . A 
comparison between the dash curve and the dash dotted curve in Fig. 2(d) shows that the concur- 
rence increases rapidly and tends to move to infinity by increasing the value of the parameter £ 2 , 
which is different from the results obtained from the Gaussian distribution and double-Gaussian 
distribution. That is to say, the strong £ 2 is helpful to keep the better entanglement for the bimodal 
distribution. 

Up to now we have examined the nearest neighbouring concurrence C(49,50) with different 
Gaussian distributions for purities and strengths of magnetic field. It is interesting to study the 
effect of the different Gaussian distributions on the concurrence for the rest of the sites in the 
chain. For the Ising model, a similar analysis can be carried out for the nearest neighbouring 
concurrence C(50,51),the concurrence is located at the centre of the double-Gaussian distribution. 
This is demonstrated in Figs. 3 and 4 by the evolutions of the concurrence. Figure 3 corresponds 
to the case in which the exchange couplings are varying, and the peak of the maximal entangle- 
ment becomes larger than that in Fig.l. It is the different distributions that lead to considerable 
different evolutions of the entanglement, hence the entanglement is rather sensitive to any small 
change in the exchange interaction for the bimodal distribution. As shown in Eq.(5), for the bi- 
modal distribution, the strengths of impurity are mostly located at two sites^g^o, 0:50,51)- The 
nearest neighbouring concurrence increases with the increasing of £ 2 , so that by adjusting ( 2 one 
can obtain a strong entanglement. The results that we have obtained here are also consistent with 
those in Ref. [20](Fig.4). Figure 4 corresponds to the case in which the external magnetic field is 
varying, the entanglement between nearest neighbours tends to be reduced in the presence of the 
external magnetic field for the double-Gaussian distribution and the bimodal distribution, while 
the entanglement between 49 and 50 increases as shown in Fig. 2. The numerical calculations also 
show that as the parameter A increases from to 4, similar behaviours to those in Figs. 2(c) and 
2(d) are shown in Figs. 4(c) and 4(d). 

From the above analysis, it is clear that the three different distributions(the Gaussian distri- 
bution, double-Gaussian distribution, and bimodal distribution) have a notable influence on the 
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nearest neighbouring concurrence. As for the case 7 7^ 1(XY model) or the next nearest neigh- 
bouring concurrence, we will present further reports in the future. 



VI. CONCLUSION 

The entanglements near particular locations in a one-dimensional s — | random XY spin sys- 
tem have been investigated. Through analyzing the exchange couplings (or external magnetic 
fields) of three different distributions(the Gaussian distribution, double-Gaussian distribution, and 
bimodal distribution), we have shown that the entanglement can be controlled and enhanced by 
varying the strengths of the magnetic field and the impurity distribution in the system. The near- 
est neighbouring concurrence exhibits some interesting phenomena. For a certain distribution, 
concurrence C(49,50) decreases with the increase of (, while concurrence C(50,51) increases. 
Different behaviours in the varying of the external magnetic field can occur close to and above the 
critical point. The different distributions play an important role in enhancing the entanglement. 
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FIG. 1 : The nearest neighbouring concurrence C(49,50) as a function of the reduced coupling constant A 
at different values of impurity ( for different distributions, with the system size N =101 and the anisotropy 
parameter 7 = 1. 
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FIG. 2: The nearest neighbouring concurrence C(49,50) as a function of the reduced coupling constant A 
different strengths of the external magnetic field £, with the system size N =101 and the anisotropy parameter 
7 = 1. 
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FIG. 3: The nearest neighbouring concurrence C(50,51) as a function of the reduced coupling constant A 
at different values of impurity ( for different distributions, with the system size N =101 and the anisotropy 
parameter 7 = 1. 
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FIG. 4: The nearest neighbouring concurrence C(50,51) as a function of the reduced coupling constant A 
different strengths of the external magnetic field £, with the system size N =101 and the anisotropy parameter 
7 = 1. 
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